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Wc investigate the dissipation of magnetic flux in primordial star-forming clouds throughout their 
collapse including the run-away collapse phase as well as the accretion phase. We solve the energy equation 
and the non-equilibrium chemical reactions in the collapsing gas, in order to obtain the radial distribution 
of the ionized fraction during the collapse. As a result, we find the ionized fraction is high enough for the 
magnetic field to couple with the gas throughout the evolution of the cloud. This result suggests that the 
jet formation from protostars as well as the activation of magneto-rotational instability in the accretion 
' (— I '' disk are enabled in the presence of the cosmological seed magnetic fiux proposed by Langer, Puget & 

PL|! Aghanim (2003). 
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INTRODUCTION 



Typical mass or initial mass function of population III stars are fundamental parameters that have great impacts 
on subsequent structure formation of the universe. Those stars are expected to be as massive as 100 — lOOOM©, which 
QQ , ionize/dissociate the surrounding media. As a result, star formation activities in the neighbourhood of the stars are 
highly regulated (e.g. Susa (2007) and the references are therein). 

There are two chief reasons that the population III stars are expected to be very massive. First one comes from 
■ numerical studies by several authors (Nakamura & Umemura 1999; Bromm, Coppi & Larson 1999; Abel, Bryan & 
f — I Norman 2000; Nakamura & Umemura 2001; Bromm, Coppi & Larson 2002). Those studies reveal that the prestellar 
. cores formed as fragments of primordial gas are as massive as ~ lO'^M© — lO^M©. Second ground is the very high 
accretion rates of those stars, which is as high as 10"'^ — lO~^M0yr~^. These two facts are the direct consequence of 
relatively high temperature (^ lOOOK) of primordial gas because of inefficient cooling by H2. 
' On the other hand, abundance pattern of the hyper metal-poor stars seems to be more consistent with that of 
JH [ "faint" supernovae as remnants of less massive (~ 25Mq) population III stars(Christlieb et al. 2004; Frebel et al. 
2005; Iwamoto et al. 2005). Recent theoretical studies suggest the possibilities for the formation of such less massive 
population III stars. If the primordial gas was once ionized, enhanced fraction of H2 causes more efficient cooling and 
HD formation. Since HD molecules have lower excitation energy than that of H2, gas can be cooled below lOOK. As 
a result, the mass of the fragments could be smaller than 100MQ(O'Shea, Abel, Whalen & Norman 2005; Johnson & 
Bromm 2006). O'Shea & Norman (2007) also demonstrates that such low mass population III stars cloud be formed 
directly from the cosmological density fluctuations. 

Those new ideas are quite interesting and promising, however, we still do not know the actual mass of population III 
stars formed from ~ IO^Mq — lO^M© prestellar cores, which are commonly found in cosmological simulations. Since 
they are only 10'* — lO^cm"'^, we have to follow the subsequent evolution of the collapsing cloud. Further evolution of 
prestellar cores is first investigated by Omukai & Nishi (1998), and they found that the collapse proceeds in a run- away 
fashion and converges to Larson-Penston type similarity solution (Larson 1969; Penston 1969; Suto & Silk 1988) with 
polytropic index 7 ~ 1.09. They also find the mass accretion rate is very large compared to the present-day forming 
stars, although spherical symmetry is assumed in their radiation hydrodynamic simulations. Recently, Yoshida et al. 
(2006) have performed cosmological simulations in which the run-away collapsing core is traced up to a very dense 
regime (nn ^ 10^"cm~^), taking the Sobolev type line transfer approximations. They also find basically consistent 
results in 3D cosmological simulations with previous ID results by Omukai & Nishi (1998). 

However, most of the mass of protostar is accumulated in the accretion phase, which has not been traced especially in 
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multi-dimensional simulations. Recent numerical simulations performed by Saigo, Matsumoto & Umemura (2005) show 
that disks would be formed at the center of collapsing primordial clouds, which might result in the disks surrounding 
the protostars. or binaries. Thus, the actual mass of a population III star should depends on the mechanism of the 
angular momentum transport in the accretion disk. There arc a few possibilities of the angular momentum transport 
mechanism such as gravitational torque by the nonaxisymmetric structures in the accretion disk, the interaction 
among the fragments of the disk (Stone et al. 2000; Bodenheimer et al. 2000), and the turbulent viscosity triggered 
by Magneto-Rotational Instability (MRI) (Hawley & Balbus 1992; Sano, Inutsuka & Miyama 1998; Sano & Inutsuka 
2001). As for MRI induced turbulence, the strength of magnetic field in the disk is the key quantity to activate 
the instability (Tan & Mckee 2004; Tan & Blackman 2004). It is also worth noting that recent MHD simulations by 
Machida et al. (2006) suggests the possibility of bipolar jets from proto-population III stars, which also could suppress 
the mass accretion onto the central core. The formation of jets also requires the presence of some level of magnetic 
field. Therefore, it is quite important to assess the magnetic field strength brought into the accretion disk from initial 
weak cosmological seed field (e.g. Langer, Puget & Aghanim (2003)). 

Maki & Susa (2004) investigated this issue by solving detailed chemical reaction rate equations coupled with energy 
equation in run-away collapsing core. They found that magnetic field is always frozen to the collapsing core in case 
the strength of initial magnetic field is smaller than 10~^(nH/10^ cm""^)"'^^ G. This is comparable to the maximal 
strength of magnetic field which allow the clouds to collapse. Therefore they conclude magnetic field is always frozen 
to the collapsing cloud in the run-away collapsing core. They also evaluate the minimal strength of magnetic field 
which is required to activate MRI assuming magnetic field is frozen to the gas not only in the run-away phase, but 
also in the accretion phase. 

This assumption is based upon the argument that the temperature of the accretion flow would rise faster than 
that of the core, because of the shock heating. As a result, the ionization degree is expected to be higher than those 
in run-away phase, which guarantee the gas to be frozen to the magnetic field lines. However, such heating is only 
important in the final phase of accretion, where the fiow settle onto the accretion disk. In order to investigate the 
coupling of matter and magnetic field in the accretion phase, we need to calculate the actual thermal and chemical 
evolution of accretion flow. Aside from the issue on ionization degree, the ambipolar diffusion velocity is proportional 
to the square inverse of density {v^^ cx p~^). As a result, the diffusion velocity increases rapidly in the mass accretion 
phase, since the density at a fixed radius decreases as the accretion proceeds. Thus, we need a detailed treatment to 
test the coupling between the gas and magnetic field in accretion phase. 

In this paper, we investigate the dissipation of magnetic field in collapsing primordial gas cloud from run-away 
phase to accretion phase, by solving detailed thermal and chemical rate equations. In the next section, we describe 
the formulations employed. In §3, results of our calculations are shown. The formation of jets and activation of MRI 
is discussed in §4. Final section is devoted to summary. 



2. METHOD OF CALCULATIONS 



In order to evaluate the coupling of the gas and the magnetic field, we need to assess the amount of ions and electrons 
in collapsing primordial cloud. The collapse of the cloud is expected to proceed in run-away fashion (run-away phase) 
in the beginning, followed by the mass accretion phase after the formation of rotationally supported disk at the center. 
In this paper, we follow the chemical and thermal evolution of the materials to form star, throughout the two phases. 
We assume spherically symmetric collapse of the progenitor gas, although we expect the formation of rotationally 
supported disk in the very dense regime (nn ^ 10^°cm~'^). It is also assumed that the magnetic field is so weak that 
the dynamics of the collapsing gas is not affected by the magnetic force. 

2.1. Run-away Collapse Phase 

The run- away collapse phase is traced by solving 1-dimensional hydrodynamics. The equation of continuity is 

dm 9 , , 

— = 4nr'p, (1) 

whereas the equation of motion is 
Du , , dp Gm 

:^_47rr2-i- — , (2) 

Dt dm 

where m is the mass within radius r, p is the density of the cloud, u is the velocity and p is the pressure. Since 
the equation of state of collapsing primordial gas in run-away phase is known to be approximated by polytrope with 
7 = 1.09 (Omukai & Nishi 1998), we employ following equation: 

p^AV, 7 = 1.09 (3) 

Wc solve above set of equations by spherically symmetric Lagrangian hydrodynamics code developed by ourselves, 
following the Piecewise Parabolic Method (PPM) described in Colella & Woodward (1984). 
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2.2. Mass Accretion Phase 

The dynamics of accretion phase could be approximated by simple free-fall, since the flow velocity is supersonic. 
Basic equations of free-falling accretion flow consists of the equation of continuity and the equation of motion: 

ATrp{tQ,ro)rldro = ATrp{t,r)r^ {t)dr, 

(fr{t) GMo 



(4) 
(5) 



df^ r2(t) ' 

where r, t and to represent the position of fluid element, time, and the time when the mass accretion starts, respectively. 
Mq denotes the mass within rg at t ~ Iq: 

Mo= / ^irr'opodr'o. (6) 

The initial radius rg is defined as tq = r{tQ). Therefore, the the solution of above equation of motion is given as 
r = r(t;to,ro). The velocity of the fluid element is given by 



GM, 







r{t) ' 



(7) 



where Sq is the total energy defined as Sq = w^(to)/2 — GMq/vq. Using equation (4), the density at {t,r) is given as 



Po 



(8) 



{dr /drojt is the partial differentiation by keeping t fixed. An explicit expression of this term is given in appendix 
1. In order to clarify the validity of the free-fall approximation, it is compared with the similarity solution with 
7 = 1.09 (Suto & Silk 1988) in accretion phase. The initial condition of the accretion flow is set as the final state of 
the run-away collapse with 7 = 1.09. Figure 1 illustrates the density and velocity distributions of several snapshots. 
It is obvious that the free-fall approximation can describe the matter distribution in accretion phase very well for 
polytropic gas. In other words, the thermal evolution has little effects on the dynamics of the accretion flow. Thus, 
the dynamics can be approximated by free-fall formula, irrespective of the internal energy equation. 





Fig. 1. Time evolution of density /velocity profiles in accretion flow are shown. Dotted lines and marks represent the result with 
free-fall approximation, whereas the solid lines denote the results from the similarity solution in accretion phase with 7 = 1.09. 



In the accretion flow, the effective polytropic index is basically unknown. Thus, we have to solve the following 
energy equation for each fluid element: 
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de d f 1 



dt ^dt\p 



C'-'"''\ (9) 



In this expression, p denotes the pressure, which is related to the internal energy e by the equation of state: 

P=(7ad-l)pe, (10) 

where 7ad is the adiabatic exponent, denotes the net energy loss rate per unit mass. 

We take into account the rovibrational line cooling by H2(>Ciinc), the continuum radiation from the gas (£cont), and 
cooling and heating associated with chemical reactions (£diss, ^h- i ^^'^ ^3body)- C^'^"^^ is the sum of all these 
contributions. 

The cooling rate by the rovibrational transition of hydrogen molecules are assessed by the fitting formula given by 
Galli & Palla (1998) in case the lines are optically thin. In the optically thick regime, we employ the escape probability 
formalism, as described in equations (12)-(15) in Omukai (2000), except that we use actual velocity gradient in our 
calculations. 

The cooling rate owing to the continuum radiation is given by 

-^^cont = 4(TKgasT'' (11) 

where a denotes the Stcfan-Boltzmann constant, Kgas is the Planck mean opacity of gas (Lenzuni, Chernoff & Salpeter 
1991), which includes bound-free absorption, free-free absorption, photodissociation, Rayleigh scattering, and collision- 
induced absorption. This formula is valid in case the accreting gas is optically thin for continuum radiation. We 
confirmed that the optical depth is smaller than unity throughout the accretion phase. 
Cooling due to the latent heat of II2 dissociation is given by 

£,,3.== 4.48^ eVs-ig-\ (12) 

where tt-h is the number density of hydrogen nuclei, and {dyn.-^ / dt) ^ is the dissociation rate of H2. 

On the other hand, the gas is heated when a hydrogen molecule is formed. The heating rate per unit mass is given 

by 

3.53 f dyii2 \ -1 



P V "Wh- l + "cr/'lH 

1.83 / dyuA ?^H _i _i 

"2 p V "W H+ 1 +"cr/nH 

4.48/d?/H2\ "H _i _i 

e^3body = —7- — -, CV S g , (15) 

P \ dt J 3,^^^y 1 + ncr/nn 

where {dyii^ / dt)ii- , {dyii^/dt)^^ and {dyii^ /dt) shody are the formation rates of H2 by process, process and 
three body reactions (HoUenbach & McKee 1979). Tier is the critical density defined as 

-1 



10" 

Ti/2 



l.eynexp 



1.4yHexp 



1.2 X lo^y 

T+ 12007 



cm"-^. (16) 



2.3. Chemical Reactions 



Since the dissipation of the magnetic flux, as well as the temperature of the gas, strongly depends on the chemical 
abundances, we have to solve the non-equilibrium chemical reaction rate equations coupled with equation of motion 
and energy equation described in previous subsections 2.1 and 2.2. 

The evolution of the fraction of species i is followed by solving the equations, 

24 24 24 24 24 

l — l m — 1 l — l m— In— 1 

(i = l,2,3,---,24), (17) 

where yt = ni/nn is the fraction of species i, kim [cw? s^^] and kimn [cm® s^^] are the reaction rate coefficients with 
respect to two-body processes and three-body processes, respectively. In our calculation, we include 24 species: c~, 
H+, H, H", H2, H+, H+, D, D+, D", HD, HD+, H2D+, He, He+, He++, HeH+, Li, Li+, Li++, , Li", LiH, 
and LiH+. We employ the latest reaction rate coefficients in the following papers, Galli & Palla (1998), Omukai 
(2000), Standi, Lepp & Dalgarno (1998), Flower (2002) and Lepp, Standi & Dalgarno (2002). As for the radiative 
recombination, we use the rate coefficients based on Spitzer (1978). These reactions are the same as our previous 
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paper(Maki & Susa 2004). Here we stress the importance to include above rare elements such as Li, because the 
coupling of magnetic field with gas can be attained by very low fractional abundance of electrons and ions (Maki & 
Susa 2004). 

2.^. Drift velocity of magnetic field 

Magnetic field is dissipated from star-forming gas via ohmic loss and ambipolar diffusion. We assess the drift 
velocity vbx of the field lines due to these two processes, which is compared to the accretion velocity of gas. We 
evaluate the drift velocity following the formulation by Nakano & Umebayashi (1986). There are two important 
quantities which characterize these diffusion processes. They are Tj/ and uji, which denote the viscous damping time of 
the relative velocity of charged particle v to the neutral particles, and the cyclotron frequency of the charged particle 

respectively. Then, t^, is expressed as 

(18) 



where /i^n is the reduced mass, n^^ rin, and p^, are, the mean number density for the charged particle the neutral 
particle n, and the mass density of charged particle v, respectively. The averaged momentum-transfer rate coefficient 
for a particle colliding with a neutral particle is expressed by {(Jv)^^ . We use the empirical formulae for the 
momentum-transfer rate coefficients (Kamaya & Nishi 2000; Sano et al. 2000). 
According to Nakano & Umebayashi (1986), the drift velocity is given by 

A c 

where 

Al 
A. 



A = Ai+Ai (20) 

2 

E Pv'''i^^v ('91 \ 

1-Kr2a'2' ' 

-'Y.Tr^.^ (22) 

^ 1 + Tf.Ujt 



B is the mean magnetic field in the primordial cloud, the suffix x means x direction component in local Cartesian 
coordinates where the z direction is taken as the direction of B. Wc replace (l/c)(j x B)x in equation (19) with the 
mean magnetic force B" /Anr, where B is the mean field strength in the cloud, r is the radius of the cloud at which 
we assess the drift velocity. 

2.5. Initial Conditions 

We consider primordial star- forming gas clouds that formed in the mini-halos with ~ 10^ Mq. We set uniform 
and spherical gas cloud with tih = lO'^ cm~^, Tq = 250 K, and M = 10* Mq. Such clouds are commonly found 
3D cosmological simulations (Bromm, Coppi & Larson 1999; Abel, Bryan & Norman 2000; Abel, Bryan & Norman 
2002; Bromm, Coppi & Larson 2002; Yoshida et al. 2003; O'Shea & Norman 2006). We also use the cosmological 
abundance given by Galli & Palla (1998), as the initial fraction of the chemical compositions. The initial magnetic 
field strength is expected to be very weak, however, its magnitude is still under discussion. Thus, we regard the field 
strength free parameter of the calculations. 

2.6. Switching from Run-away Phase to Accretion Phase 

The run-away phase and the subsequent accretion phase have to be treated in different formulation from each other 
as described in subsections 2.1 and 2.2. Thus, we have to switch the scheme from the method in run-away phase to 
that in accretion phase. In the light of physical arguments, two phases should be switched when the central run-away 
collapse is stopped due to thermal pressure or centrifugal force. In our calculation, we switch the scheme when the 
sonic point of the accreting flow go inside a certain small radius. Wc set this radius to be fsw = 700 Rq, which is 
comparable to the disk size when the accreted mass is comparable to ^ 1 A/0 (Tan & Mckee 2004). It is also worth 
noting that the central density of the core at switching is 



3. RESULTS 



3.1. Runaway Collapse Phase 

3.1.1. The distribution of the chemical abundances 

Two snapshots of yi distribution during the run-away collapse phase are shown in Figure 2. The left panel shows 
the distribution at the time when the central density tzh.c reaches 10^^ cm""^, whereas the right panel represents the 



6 



Maki & Susa 




10' 
10 
10' 
10" 
10"' 
10"' 
10' 

io" 

10 

10 

10 

10 
10 
10 
10 
10 
10 
10 



-10 



-15 



I I I I I I 

- lo'^ [cm ^ ] 
H ^ 




•^2 

iiiNi I 



iiii / 



10 10 10 

r [cm] 



io" lo" io'^ lo" io" io" 

r [cm] 



Fig. 2. The distribution of the main species, e , H+, H, H2, HJ, H^, Li, and Li+ are plotted. The vertical axis denotes the 
fractional abundance yi of the above species, and the horizontal axis denotes the radius from the core center. Left panel represents 



the snapshot when the central density satisfies ?ih,c = 10 cm , whereas right panel shows the results at ?ih,( 
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Fig. 3. Drift velocity vbx as a function of the radius r and the field strength B at the central density . Contour map of 
v^^/u is shown on r — B plane for nH,c = lO'^^ cm~^ Solid curves represent the constant loci along which log{vg^/u) equal 
to the values labeled on the curve. The dashed curve represents the critical field strength Bcr given by the equation (23). 



distribution at tt-h.c = 10^^ cni""^. The horizontal axis denotes the radius from the center of the primordial gas cloud 
r [cm], and the vertical axis is the fractional abundances yi of each species. Since the central density increases as the 
collapse proceeds, right panel represents the later epoch than that of the left. 

Roughly speaking, in both of the panels, fraction of electrons (ye) in the inner region are smaller than that in 
the outer part. These results could be simply understood since the recombination process in the inner dense region 
proceeds faster than in the outer envelope. We also find the opposite behavior in the inner most region of right panel. 
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at which j/o increases as r decreases, because of the coUisional ionization. It is also worth noting that Li becomes 
the main provider of electrons around r ^ lO^^'cm at ?iH,c = 10^® cm^"^. Because of the combined effects of these two 
(collisional ionization & presence of Li), yc never gets lower than 10~^^ as far as we consider njj.c < lO""^^ cm"'^ in 
run-away phase. 

3.1.2. Drift velocities in run-away collapse phase 

In Figure 3, we show the ratio of the drift velocity v-qx to the infall velocity u when nH,c = 10^^ cm~'^ is satisfied. 

The vertical axis shows the magnetic field strength, whereas the horizontal axis shows the distance from the cloud 
center. The solid curves in Figure 3 show the contours of \ogiQ{vBx / u) . The magnetic fields are dissipated in the 
region where \ogiQ(vBxlu) > 0, in contrast, the fields are frozen to the gas in the region where log]^Q(?;Ba:/u) < 0. The 
dashed curve in Figure 3 is the critical field strength B„ that is defined by the equation 

Bi^GMirprl^ (23) 
Airr r^ 

Note that since we are interested in the collapsing gas cloud, the magnetic force needs to be weaker than the gravi- 
tational force. Hence, our calculations are valid in the region where the field strength in the cloud satisfies B < Ba-. 
We find clearly from Figure 3 that the frozen-in condition vbx/u < 1 is almost always satisfied as long as B is less 
than the critical field strength Bcr- We also find the basically same results for other snapshots, i.e. the drift velocity 
is smaller than the infall velocity anytime and anywhere if _B < Ser- 
in addition, if i? < Bct is satisfied at some initial time and position (to,fo); this condition also holds at some later 
epoch {t,r{t)). This statement is proved as follows: Combining the magnetic fiux conservation equation 

2TrrodroB{to,ro) = 27Tr{t)drB{t,r{t)) (24) 

and the mass conservation law given in equation (4), we obtain 

i?(t,r(0)=i?(^o,.o)^%:^, (25) 
Using equations (23) and (25) we have 



B{t,r{t)) _ S(io,ro) f p{t,r{t)y (t) 



B„{t,r{t)) B„{to,ro) \ pito,ra)rl 



1/2 



(26) 



The second term in right side is always less than or equals to unity in the run-away collapsing cloud, since p oc r^^ 
in the envelope, and p (x in the core. Thus, if B{tQ,rQ)/ Bci{to,ro) < 1 is satisfied, B is less than _Bcr throughout 
the collapse. Remark that we assume the fiux conservation as equation (25), which gives the maximal field strength. 
Therefore, above inequality {B < Ba-) also holds even if the magnetic flux is dissipated. 

In summary, it is concluded that the magnetic fields are always frozen to the whole cloud in the course of run-away 
collapse phase, if the magnetic force is much smaller than gravitational force at the beginning of the collapse. 

We also remark that the slight increase due to collisional ionization at nn.c = 10^^ cm~'^ implies that j/e will become 
larger as the central density increases up to nn.c > 10^^ cm~^, since the temperature gets higher as the density increases 
for 7 ~ 1.09 polytropic gas. As a result, the magnetic field is expected to be frozen to the gas even for nn.c > 10^^ cm"-^ 
in accretion phase. 

3.2. Accretion Phase 

3.2.1. Evolution oj gas and chemical species in accretion phase 

The evolutionary sequences of number density nn, temperature T, velocity u, and electron fraction y^ in the accretion 
phase are illustrated in Figure 4. 

Upper left panel illustrates the evolution of density profile. The solid curves labeled as 0-4 corresponds to the epoch 
at which the central accreted mass equals to 8.87 x 10~^ M0,1.58 M0,12.7 M0,5O.2 Mo,100 M©, respectively. It 
is clear that the density for r ^ decreases as the collapse proceeds, which is basically the same feature found in 
similarity solution of polytropic cloud in accretion phase. 

Because of the high accretion velocity at later snapshots (lower left), the gas is heated efficiently by adiabatic 
compression. As a result, temperature gets higher as the collapse proceeds, except at the beginning of the accretion 
phase, when the gas is cooled efficiently due to the enhanced fraction of II2 molecules. 

The ionization degree (electron fraction y^, lower right) also increases as the collapse proceeds (Figures. 4 and 5), 
due to the higher temperature for later snapshots. As a result, y^ never gets lower than 10~^^. We also find that Li 
is not so important as was in the final phase of run-away phase, since the hydrogen win back the position of chief 
provider of electrons in the accretion phase( Figure 5). 
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Fig. 4. The evolutionary sequences of the primordial gas cloud in accretion phase are shown. Four panels (a)-(d) show the spatial 
distribution of following physical variables as functions of radius: (a) number density nji, (b) temperature T, (c) velocity u, and 
(d) electron fraction j/c- Five time sequences (0-4) are plotted. Accreted central mass Mc{t) is used as a clock. Corresponding 
mass at the stages are: 0:A/c(i) = 8.87 X 10"^ Mq, l:Mc{t) = 1.58 Mq, 2:Mc{t) = 12.7 Mq, 3:Mc{t) = 50.2 Mq, 4:Mc(t) = 100 Mq. 



3.2.2. Drift velocities in accretion phase 

As briefly discussed in section 1, the ambipolar diffusion velocity is proportional to the inverse square of gas density. 
Its dependence on various physical quantities are described as oc y~^p~^B^r~^ (see equations (19)-(22)). On the 
other hand, the accretion velocity scales as u oc {Mjry/^ for r ^ 0. Since we consider the accretion phase starting 
from the final phase of run-away collapse with 7 = 1.09, the density profile and the central mass in accretion phase 
also depend on 7. According to Suto & Silk (1988), we have 

p(i, r) oc t(2-37)/2^-3/2 for r ^ (27) 
Mc(t)oct4-37 (28) 

Thus, the limiting behaviour of the infall velocity is described as 

Uit, r) oc i(4-37)/2^-l/2 foj. ^ ^ Q (29) 

Combining above set of equations, we have the dependence of the ratio v^^^ /u on t: 

amb 

^ OC y-ls2t(97-8)/2^5/2 (30) 

This equation indicates that the ratio keeps growing in accretion phase for fixed r,B and yc, since 97 — 8 is positive 
for 7 = 1.09. In reality, however, j/c increases rapidly as the accretion proceeds (see Figures. 4 and 5), and it offsets the 
increase of v^^^ /u. 

Figures 6 and 7 illustrate the contours of logrg(?;Bx/w) on r — B plane for two epochs. Two figures correspond to 
the snapshots when the central accreted mass satisfies = 12.7 Mq and M^. = 100 Mq, respectively. The notations 
are same as Figure 3. It is clear that the drift velocity is always smaller than the accretion velocity in two snapshots. 
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Fig. 5. The radial distributions of the fractional abundances of the main species, e, H+, H, H2, H , Hg , Li, and Li+ 
in the accretion phase. Two panels correspond to the snapshots at Mc(t) = 12.7 M0(left), and Mc(t) = 100 MQ(right). 
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Fig. 6. The same as Fig. 3, except that the contours arc plotted in the accretion phase characterized at Mc =12.7 Mq. 

as long as B < Ba- In fact, we find that this is true all through the calculations. Besides, the equation (26) holds also 
in the accretion phase. Considering that the density profile in the accretion phase is poc r~^-^ ^, B is less than Bcr 
throughout the collapse, if B{to,ro)/ Bci{to,rQ) < 1, as discussed in §3.1.2. In other words, magnetic force is always 
negligible if it can be ignored in the beginning of the collapse. 

Thus, it is concluded that the dissipation of magnetic flux is negligible throughout the mass accretion phase, as well 
as the run-away collapse phase. 



4. DISCUSSION 



We confirmed that the magnetic field is frozen to the star-forming primordial gas cloud even in the accretion phase. 
Here we discuss the possibility of jet formation and activation of MRI considering the magnetic field strength brought 
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Fig. 7. The same as Fig.6, except Mc = 100 Mq. 



into the accretion disk surrounding the protostar. 

The magnetic field strength brought into the accretion disk is assessed under the frozen-in condition as follows: 

2 



R 

'"disk 



(31) 



Here Bq is the initial field strength, rjisk denotes the disk radius, whereas R describes the initial radius within which 
includes the total mass of the disk-star star system M^:^is'k'- 

i^=(3A4disk/4^Po)'/^ (32) 
where po represents the initial density of the cloud, rdisk can be evaluated by equation (15) in Tan & Mckee (2004): 



2 / A/- \^/'^ 
-^^'-'+disk 



(33) 



rdisk ^ 66.4AU f ^\ , 

where /Kep denotes the ratio of the rotation velocity to the Kepler velocity of accreting matter, which is found to be 
^ 0.5 in Abel, Bryan & Norman (2002). Thus, we have the magnetic field strength in the disk as follows: 



B, 



disk - 



:7.5x 10^i"G 



So 



3.7x 10-16G 



( \ 
Vl03cm-3y 



0.5 



disk 



-40/21 



(34) 



lOAfo . 

Several possibilities to generate cosmological seed magnetic field have been proposed so far. Most of the mechanisms 
predict i?iGM ^ 10~^^G(Widrow 2002), except the magnetic field generated by radiation transfer efi'ects of powerful 
ionizing sources such as quasars or first stars (Langer, Puget & Aghanim 2003). They suggests the possibility to 
generate coherent magnetic field with Sigm ~ 10~^^G. Since the magnetic field is frozen to the primordial gas at low 
densities (tt-h < lO'^cm"^), Bq at nn = lO^cm"^ can be evaluated as 



:3.7x lO^i^G 



BlGM 

10-19G ; Vl03cm 



2/3 / I + z 

"20" 



(35) 



Recently, 3-dimensional MHD simulations on primordial star formation have been performed by Machida et al. 
(2006), assuming ideal MHD condition is always satisfied in the collapsing gas. In fact, our present results guarantee 
this hypothesis. They found that the protostellar jet is driven in primordial environment if Bq ^ 10~^G at nn = 
lO'^cm"'^. Comparing this condition with equation (35), it is concluded that 1) jets could be driven in first star 
forming clouds if the seed field is generated by the mechanism proposed by Langer, Puget & Aghanim (2003), 2) 
whereas the other mechanisms cannot generate the seed field enough to drive the jets. 
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MRI can be activated in the accretion disk in case the magnetic field in the disk is larger than a critical value (Tan 
& Blackman 2004): 

This threshold is assessed by the confrontation between the growth rate of MRI and the ohmic dissipation rate. 
Combining equations (34)- (36), we find that MRI is activated in case the seed field satisfy 

-2 / J \ 4 / , r \ 2.155 



B 



IGM 



>,5x:o-"G (i^)"(f-) (^)- 

-3/4 . X 1/2 / , xl/2. ^ x-3/4 



T 



(37) 



^lO^Ky V 10 y V5 X 10-1" gcm-'V V6OOi?0, 

Therefore, MRI is driven only if the seed field generation mechanism by the transfer effects of ionizing radiation works. 
Based upon these arguments, we emphasize that the mechanism proposed by Langer, Puget & Aghanim (2003) should 
be scrutinized since their results still based upon the argument of order-estimation. 



5. SUMMARY 



In this paper, we investigate the dissipation of magnetic flux in star-forming primordial gas cloud. We solve non- 
equilibrium chemical reaction equations, coupled with thermal and dynamical evolution of the collapsing cloud all 
through the run-away phase as well as the mass accretion phase. Thus, we obtain the detailed evolution of ionized 
fraction of the gas, which enables us to assess the coupling between gas and magnetic field. 

As a result, we find that the magnetic field is basically frozen to the gas anywhere in collapsing star-forming 
primordial clouds at any time. Based upon this result, we find the cosmological seed magnetic field generated by most 
of the mechanisms proposed so far is not sufficient to form jets as well as to activate MRI in the star-forming cloud. 
Only one mechanism proposed by Langer, Puget & Aghanim (2003) is able to create sufficient field strength. 

We thank Kazu Omukai for stimulating discussions. Noriaki Shibazaki and Ken Ohsuga are acknowledged for 
continuous encouragement. The analysis has been made with computational facilities at Rikkyo University. This work 
was supported in part by Ministry of Education, Culture, Sports, Science, and Technology (MEXT), Grants- in- Aid, 
Specially Promoted Research 16002003 and Yoimg Scientists (B) 17740110. 



Appendix 1. Density distribution of free- falling matter 



/ dr 

In this Appendix, we derive the explicit formula for — — 1 which gives the density distribution of free-falling 



matter in equation (8). 

Equation of motion (5) is integrated as 

1.2 GMq GMo 1 2 ^ n^ /Ai\ 

7;r = to, to = 77^0, (to > 0), (Al) 

2 r To 2 

where uq is the velocity of a fluid clement at (to,''o) where tq denotes the position of the element at some initial time 
tQ. A solution of this equation is given as 

r{a,ao)^ ^-^{l-cosa), (A2) 
Zto 

GAI 

t{a, ao) = to- [{a - sin a) - {ao - sinag)] , (A3) 

where a is the so-called development angle, and ao is its value at i = to- ao also satisfies following relation 

GMq s / a ^ \ 

-Z7^[i-cosao). (A4) 
Using above relations, (9r/9ro)t is derived as 
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